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5 1. Claim of Priority. 

This patent application claims priority to European Patent Application serial number 
02021903.6 filed on September 27, 2002. 

2. Field of the Invention. 

This invention relates to the field of signal processing, and in particular to a parametric 
1 0 recursive digital filter. 

3. Related Art. 

As known, digital filters include at least one delay unit and a coefficient network, which 
determine the cut-ofi7center frequency (i.e., f c ) of the filter. In some application it is desirable to 
change the cut-off or centre frequency of a filter in response to a control quantity (parametric 
15 filter). 

In conventional digital filters, the cut-off or center frequency is controlled with either the 
aid of coefficient look-up tables or real-time calculation of the filter coefficients. However, 
known structures for implementing the coefficient look-up tables and/or the real-time calculation 
require a rather large amount of memory and/or increased computational ability. This is of 
20 course particularly true in the case of higher-order filters. W. Schussler, W. Winkelnkemper, 
suggest in the publication entitled "Variable Digital Filters", Arch. Elektr. Ubertr., Vol.24, 
1970, issue 11, pages 524-525, an arrangement in which the delay elements of a conventional 
digital filter are replaced by all-pass filters in order to obtain a parametric filter. However, the 



authors of this publication concede that this procedure cannot be implemented in digital filters 
because the all-pass filters comprise an attenuation when the highest all-pass filter coefficient is 
not equal to zero. As a solution for this problem, it is suggested either to use only all-pass filters 
in which the highest all-pass filter coefficient is equal to zero or only to use non-recursive filter 
5 structures. A problem with the proposed solution is that few filters can be implemented so that 
the possible applications are very limited. 

Therefore, there is a need for a parametric recursive digital filter. 

Summary 

10 The cut-off frequency or center frequency of a digital filter may be controlled by 

changing a coefficient value within a delay unit of the filter. A linear relationship exists between 
the coefficient value and the cut-off or center frequency of the recursive digital filter. 

The digital filter may be implemented with a recursive filter structure, such as for 
example digital wave filter structures. A filter element with controllable phase angle (e.g., 

1 5 all-pass filter ) is provided as the delay unit. 

In complex filters and higher-order filters, a plurality of delay units may be used, and the 
delay units formed by identical filter elements selected in the same manner. As a result, the filter 
elements have the same coefficient value so coefficient calculation only needs to be performed 
once and can be used for all the filter elements. 

2 0 In one embodiment, the delay unit may include an all-pass filter. The delay unit may 

include a first adder, one input of which forms the input of the delay unit, and a second adder, the 
output of which forms the output of the delay unit. A coefficient section is connected between 
the output of the first adder and a first input of the second adder. A first delay element is 



2 



connected between the input of the delay unit and a second input of the second adder, and a 
second delay element is connected between the output of the delay unit and a second input of the 
first adder. The phase angle of the filter element may be adjusted by changing the coefficient of 
the coefficient section, and the output of the first and/or second delay element being provided for 
5 connecting a feedback path. In addition to this all-pass filter embodiment, other embodiments of 
an all-pass filter may also be used. 

To simplify the structure of one embodiment of the all-pass filter, two delay units 
comprising delay elements can be interconnected with one another (in each case) in such a 
manner that only a total of three delay elements are provided, one delay element being used for 
1 0 both delay units. 

Other systems, methods, features and advantages of the invention will be, or will become, 
apparent to one with skill in the art upon examination of the following figures and detailed 
description. It is intended that all such additional systems, methods, features and advantages be 
included within this description, be within the scope of the invention, and be protected by the 
1 5 following claims. 

Description of the Drawings 

The invention can be better understood with reference to the following drawings and 
description. The components in the figures are not necessarily to scale, emphasis instead being 
2 0 placed upon illustrating the principles of the invention. Moreover, in the figures, like reference 
numerals designate corresponding parts throughout the different views. 

Figs. 1 A and IB are block diagram illustrations of two first-order all-pass filters; 
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FIG. 2 is a graph illustrating the phase response/frequency bending function, of a first- 
order all-pass filter as a function of the all-pass filter coefficient; 

FIG. 3 is a graph illustrating the frequency resolution of a parametric filter for various all- 
pass filter parameters y, 

5 FIG. 4 is a block diagram illustration of a prior art finite impulse response (FIR) filter; 

FIG. 5 is a block diagram illustration of a parametric finite impulse response filter; 
FIG! 6 is a detailed block diagram illustration of a computationally efficient structure for 
the parametric FIR filter illustrated in FIG. 5; 

FIG. 7 is a block diagram illustration of a prior art biquad filter as an example of a 
1 0 second-order infinite impulse response (IIR) filter; 

FIG. 8 is a block diagram illustration of the structure of a biquad filter as an example of a 
second-order infinite impulse response (IIR) filter; 

FIG. 9 is a block diagram illustration a parametric infinite impulse response filter with a 
low-pass filter as the first delay unit; 
15 Figs. 1 OA- IOC illustrate a series of coordinate patterns illustrating the pole-zero shift in 

the Z domain as a function of the all-pass filter parameter 7 in a second-order parametric IIR 
filter structure; 

FIG. 1 1 is a graph illustrating the pole-zero diagram of a normalized filter; 

FIG. 12 is a graph illustrating the amplitude/frequency response of the normalized filter; 
2 0 FIG. 13 is a detailed block diagram illustration of a computationally efficient 

embodiment of the second-order parametric IIR filter illustrated in FIG. 9; 

FIG. 14 is a graph illustrating the impulse response of a target filter for a second-order 
parametric IIR filter; 
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FIG. 15 is a graph illustrating the resultant amplitude/frequency response of the target 
filter; and 

FIG. 16 is a graph illustrating the pole-zero diagram of the target filter. 

5 Detailed Description 

A digital filter includes delay elements configured and arranged as all-pass filters (e.g., 
first-order low pass filters). The all-pass filter may be referred to as a frequency distorting filter 
since the specifications of the filter are often implemented on a distorted non-uniform frequency 
axis. 

10 A shift in the amplitude/frequency response |H(z)| of the filter can be achieved if z is 

replaced or mapped by an expression W(z), where W(z) meets the following conditions: 

(i) the inside of the unit circle (Z plane) must be mapped again to the inside of the 
inner circle; and 

(ii) the unit circle must be mapped to itself. 

15 

Condition (i) specifies that if H(z) is stable, H(W(z)) is also stable. The second condition makes 
it possible to shift (i.e., "map") the amplitude/frequency response into another arbitrary 
frequency range. 

W(z) is the transfer function for an all-pass filter. The amplitude/frequency response 
2 0 |H(z)| of the all-pass filter is constantly equal to one, which ensures that the amplitude/frequency 
response of the overall filter is not changed and that the first condition (i) is satisfied. Thus, a 
displacement or bending of the frequency axis occurs, and the displacement is controlled by the 
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filter coefficient of the first-order all-pass filter. As a result, filter coefficient value y of the first- 
order all-pass filter is also called frequency curvature parameter. 

FIGs. 1A and IB are block diagram illustrations of first-order all-pass filters. FIG. 1A 
illustrates a first order all-pass filter 20 that receives an input signal sequence x[n] on a line 22. 
A summer 24 receives the input signal on the line 22 and the past value of the output signal y[n- 
1] on line 26, and provides the resultant sum on line 28. The sum is input to a coefficient section 
30, which multiplies the summed signal on the line 28 by coefficient value 7. The product is 
output on line 32 to summer 34, which also receives a signal on line 36 indicative of the past 
value of the signal on the line 28. The summer 34 provides output signal y[n] on line 38. 

FIG. IB illustrates another first order low pass filter 40. This filter receives the input 
signal x[n] on line 42 and sums this signal with a coefficient weighted feedback signal on line 
44. A summer 46 provides a summed value on line 48, which is input to a coefficient multiplier 
50 that multiplies the signal on the line 48 with coefficient value -y, and a delay element 51. 
The resultant product is output on line 52 and summed with a delayed version of the signal on 
the line 48, to provide output signal y[n] on line 54. 

A first-order all-pass filter has the following transfer function: 



The first-order all-pass filter changes its phase response as a function of the all-pass parameter y, 
which can move within a range of -Ky<+1. The change of this frequency response 
corresponds to bending of the frequency axis or mapping into a new frequency range. 



FIG. 2 shows the frequency response/frequency bending function of the first-order all- 
pass filter as a function of the all-pass coefficient/frequency curvature parameter y. Phase in 
degrees is plotted along the vertical axis, while normalized frequency (i.e., frequency with 
respect to the sampling frequency f s divided by two) is plotted along the horizontal axis. As 
5 shown, high frequencies are mapped to low frequencies for negative y values, and low 
frequencies are mapped to high frequencies for positive y values. The greater the amount of y, 
the greater the frequency bending. There is no frequency bending for the coefficient value y = 0 
- that is, the frequency axis is mapped to itself again in this case (bisector). 

Since delay elements in traditional filter structures do not appear in the form of z but in 
10 the form of z \ it is also appropriate for filters with frequency bending to work not with W(z), 
but with W _1 (z). In this way, it is impossible to replace all delay elements z" 1 with their 
frequency-dependent counterpart: 



W'(z) = D(z) = -^ r EQ.2. 

15 

Due to the bending of the frequency axis, some frequency ranges are resolved more finely 
and others more coarsely in dependence on the frequency curvature parameter y. If a filter has a 
local frequency resolution of Af(f) before the "bending", the bent filter has a local frequency 
resolution of Af w (f), where the following holds true: 

20 

a/.(/) = '' C,.^, E( 3- 3 

1 + y 2 + 2 * y * cos( ) 

Js 
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where f s is the sampling frequency and f w (f, y) stands for a function that specifies the 
corresponding frequency in the bent frequency range for a given frequency and all-pass 
parameter. The frequency of the parametric filter as a function frequency /and the all-pass 
5 parameter y can be expressed as: 



. . .2V/, 
X*sm( — ) 

fM,r) = f+-* arctan ( EQ - 4 

n \- y *cos( - J ) 

Js 



The turning point frequency ftp is the frequency that the frequency resolution of the parametric 
1 o filter is the same as before the frequency bending. It can be shown that the following holds true 
for ftp: 



f = -^-*arccos00 E Q- 5 



15 For frequencies less than ftp, the frequency resolution is improved with a positive y. The 
frequency resolution decreases for negative y. FIG. 3 illustrates the frequency resolution of the 
parametric filter for various all-pass parameters y. The turning point frequency ftp is given when 

the frequency resolution is equal to one. 

FIG. 4 is a block diagram illustration of a prior art finite impulse response (FIR) filter. 
20 The filter includes delay elements 62, 64. A frequency-bending FIR filter is obtained by 
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replacing the delay elements 62, 64 (FIG. 4) with frequency-dependent delay elements 66, 68 
having a transfer function D(z), as shown in FIG. 5. 

FIG. 5 is a block diagram illustration of a parametric finite impulse response (FIR) filter 
65. If the frequency-dependent delay sections 66, 68 are replaced by their corresponding all-pass 
5 filter structure (e.g., by the structure shown in FIG. 1A), and the structure is simplified by 
dissolving redundant branches, a computationally efficient structure of a frequency-bending FIR 
filter is obtained as illustrated in FIG. 6. 

FIG. 6 is a detailed block diagram illustration of a computationally efficient embodiment 
for the parametric FIR filter illustrated in FIG. 5. The filter includes a plurality of series- 
10 connected delay elements (e.g., 70-72), summing nodes 74,75 each inserted between two delay 
elements, and coefficient sections 78, 80. The input signals for the coefficient sections 78, 80 are 
picked up at the input of the delay element 70, the output of the summing node 74, and output of 
delays 71, 72. The coefficient sections 78, 80 include programmable/adjustable coefficient value 
7. The coefficient sections 78, 80 provide coefficient section output signals on line 82, 84, 
15 respectively, which are summed with associated delay element output signals to provide output 
signals to a multiplication network 90. Products from the multiplication network 90 are input to 
a summer 92, which provides an output signal on line 94. 

Replacing a delay element in one of the traditional IIR filter structures with a frequency- 
dependent (dispersive) delay element as has already been explained for the FIR filter provides a 
2 0 structure that contains feedback branches are free of delay, but which cannot be implemented in 
this way, as has been explained in the prior art mentioned initially. FIGs. 7 and 8 show a 
recursive biquad filter before and after the mapping, respectively. FIG. 7 is a block diagram 
illustration of a prior art biquad filter as an example of a second-order infinite impulse response 



9 



(IIR) filter. FIG. 8 is a block diagram illustration of a modified biquad filter as an example of a 
second-order infinite impulse response (IIR) filter. The structure of a frequency-bending IIR 
filter (i.e., frequency warping IIR filter - WDR) is obtained by replacing the delay elements 
illustrated in FIG. 7 with frequency-dependent delay elements 80, 82 having a transfer function 
D(z), as shown in FIG. 8. FIG. 7 is based on the arrangement illustrated in FIG. 4, and has been 
expanded by a feedback network. 

FIG. 8 illustrates a WIIR filter that is suitable for higher order filters, and provides for all 
other delay elements apart from the first one to be replaced by first-order all-pass filters D(z). 
The first delay is replaced by a first-order low-pass filter H(z). This measure makes it necessary 
to recalculate (to map) the filter coefficients of the original IIR filter to the new structure, which 
is done by a relatively simple recursive formula. FIG. 9 shows the corresponding frequency 
mapping/warping IIR filter structure that includes a low-pass filter H(z) 100 and a plurality of 
first-order all-pass filters D(z) 102, 104. 

The digital filters according to an aspect of the present invention include delay elements 
configured and arranged as all-pass filters (e.g., particularly of the first order). However, to 
avoid the zero-delay feedbacks, feedback is not effected after the respective dispersive delay 
units (in the prior art) but, according to the invention, directly after the delay element within a 
dispersive delay unit. The forward branch (FIR part of the filter) with its coefficients b 0 , bi, b 2 
etc. remains unaffected by this measure. We shall now discuss a procedure for designing a 
second-order frequency-bending IIR filter with the design of a 10 Hz high-pass filter. 

Using frequency-bending filters, a fixed "normalized filter" can be designed that has an 
arbitrary fixed cut-off frequency f c , and which can be displaced into any desired frequency range 
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with the aid of the frequency curvature parameter y. FIGs. 10A-10C illustrate the pole-zero shift 
in the Z domain as a function of the frequency curvature parameter y. 

Since originally filters are assumed that are designed in the direct form, it is appropriate 
to utilize the advantages of the direct form as far as possible. The frequency range within which 
the direct form has particularly good characteristics in the sense of quantization sensitivity and 
dynamic range, is around f s /4. For this reason, the cut-off frequency f c of the prototype filter is 
left at f c = f s /4, where f s is the sampling frequency. FIGs. 1 1 and 12 show the associated pole- 
zero diagram and the amplitude/frequency response of the normalized filter, respectively. 

The frequency curvature parameter y is be calculated so the normalized filter is shifted 
towards the required cut-off frequency. For this purpose, the equation: 



2*n* f 
y * sin( ) 

f P M,r) = / + ~ * ar*an( ^7") E ^ 6 

I- y* cosier- 1 -) 
Js 



can be transformed in such a manner that the frequency curvature parameter y can be calculated 
from it. The calculation formula that produces the frequency curvature parameter y with the aid 
of the fixed cut-off frequency of the normalized filter f pre and the desired cut-off frequency f c can 
be expressed as: 



r = -( 



tan( 

J s 



) 



sin( -^-) + cos( — ^) * tan( ) 

fs Js J* 



EQ.7 
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If f pre is left at f s /4, the previous equation is reduced to: 

y^^LlL) EQ.8 
4 f s 

The expression f pre for the fixed cut-off frequency of the normalized filter was selected since it 
can be said that this filter is designed in the distorted frequency range. One can return from the 
distorted frequency range into the original frequency range by operating the parametric filter 
with the frequency curvature parameters y calculated above. 

FIG. 13 illustrates a computationally efficient embodiment of a second-order parametric 
IIR filter. This embodiment is obtained by expanding the arrangement according to FIG. 6 by 
one feedback network. From the second-order parametric IIR filter, the impulse response shown 
in FIG. 14 and the resultant amplitude/frequency response of the target filter shown in FIG. 15 is 
obtained. The associated pole-zero diagram of the target filter is shown in FIG. 16. 

The digital filters may be implemented in microprocessors, signal processors, 
microcontrollers, computing devices et cetera. The individual filter components such as, for 
example, delay units, delay elements, coefficient sections are then hardware components of the 
microprocessors, signal processors, microcontrollers, computing devices, et cetera, that are 
correspondingly used by the executable software. 

The illustrations have been discussed with reference to functional blocks identified as 
modules and components that are not intended to represent discrete structures and may be 
combined or further sub-divided. In addition, while various embodiments of the invention have 
been described, it will be apparent to those of ordinary skill in the art that other embodiments and 
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implementations are possible that are within the scope of this invention. Accordingly, the 
invention is not restricted except in light of the attached claims and their equivalents. 
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